Magnetism of Cold Fermionic Atoms on p-Band of an Optical Lattice 
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We carry out ab initio study of ground state phase diagram of spin- 1/2 cold fermionic atoms 
within two-fold degenerate p-band of an anisotropic optical lattice. Using the Gutzwiller variational 
approach, we show that a robust FM phase exists for a vast range of band fillings and interacting 
strengths. The ground state crosses over from spin density wave state to spin-1 Neel state at half 
filling. Additional harmonic trap will induce spatial separation of varies phases. We also discuss 
several relevant observable consequences and detection methods. Experimental test of the results 
reported here may shed some light on the long-standing issue of itinerant ferromagnetism. 
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Introduction: Recently, research of ultracold atomic 
gases have simulated a new wave of studying the many- 
body problems. One can create periodic potentials to 
confine the ultracold atoms by intersecting laser beams. 
Because of the experimental controllability of atom num- 
ber, dimensionality, geometry as well as interaction 
strength etc., the optical lattice provides an ideal play- 
ground for investigating many body problems. Specif- 
ically, the experimental realization of a Mott-insulator 
phase transition in the optical lattice [H0 brings us to 
the forefront of the strongly correlated systems. More 
recently, the experimental progress makes it possible to 
put fermions into the optical lattices and control the in- 
teraction between them These advances open a 
new channel to investigate numerous phenomena that 
play important roles in the condensed matter physics, 
such as quantum magnetism [f| [|| [3] , high-temperature 
superconductivity[8| and physic associated with band de- 
generacy [9(1 [ToT] . 

Itinerant ferromagnetism in transition metals is an- 
other controversial issue in condensed matter physics. 
Despite long history dating back to Stoner 1 121 . a fully 
understanding of its mechanism is still lacking [l3j . How- 
ever, a consensus is reached that the appearance of robust 
ferromagnetism is not a generic feature of the single band 
Hubbard model in a cubic lattice. Non-bipartite lattice 
structures may stabilize FM order, degenerate bands and 
the Hund' rule coupling may also help. Testing various 
scenarios using the ultra-cold fermionic atoms in optical 
lattice may shed light on this issue. The highly tun- 
able and neatness of cold atom experiments also permit 
calculations with no adjustable parameters, the phase 
diagram and experimental signatures could be deter- 
mined through the following first principle parameters: 
the depth of the optical lattice V for each directions and 
the s-wave scattering length a s . 

In this paper, we first perform band structure calcula- 
tions to identify a parameter region of V and a s where 
fermionic atoms on an optical lattice could be effectively 
described by a two-band Hubbard model. Band calcu- 
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FIG. 1: (a). Band structures of four low lying bands (s, p x , 
p y and p z ) of an optical lattice with V x = V y = 16Er, V z = 
36Er. (b)-(d). Projected density of states (PDOS) for each 
band. The bandwidth of p x and p y band (c) is about 10 times 
wider than the s-band (b), due to the large overlap of p-band 
atomic wave function. The low lying s-band shows 2D char- 
acterizes, while PDOS of p x and p y band are degenerate and 
show quasi- 1 dimension characterizes. A weakly breaking of 
the particle-hole symmetry also could be recognized from the 
asymmetry of the PDOS. Interestingly, the PDOS of p z -band 
has 3D characterizes (d), this is due to the strength of cr-bond 
along z-axis is strongly suppressed by large V z and becomes 
comparable to the weak 7r-bond along x and y directions. 



lation also provides effective coupling strengthes of the 
model. We then obtain zero temperature phase diagram 
as function of the coupling strength and the density of 
fermionic atoms by applying Gutzwiller variational ap- 
proach. Finally, we consider the effect of external har- 
monic trap within the local spin density approximation 
(LSDA) , which is shown to induce concentric shell struc- 
tures of different phases. In particular, a robust fer- 
romagnetic shell composed of spin polarized fermions 
emerges. 

The model and the band structure calculations: Six 
counter-propagating lasers generate potential of the form 
V(r) = V x sin 2 (k L x)+V y sin 2 (k L y)+V z sin 2 (k L z). Atoms 
trapped by this potential form a simple cubic lattice with 
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lattice spacing a 



Fig Q] shows the band structures 



of the optical latices with V x — V y = X&Er, V z = 36-Er, 



where Er = 



h-ki 



is the recoil energy. The degenerate 



p x and Pj,-band are well separated from the low lying 
s-band, p z band are pushed up by the large V z . As- 
suming the s-band is fully occupied, as long as the band 
gaps are much larger than the interaction scales within 
each band, we can focus on the degenerate p x and p y 
bands. Authors of 14] assume that atoms mainly inter- 
act in the p-wave channel thus could be effectively treat 
as spinless fermions and they studied effects of differ- 
ent lattice structures. In the present study, we focus on 
cubic lattice and consider the s-wave pseudo-potential 
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S(r — r'), which makes the problem 



more relevant to condensed matter system where spin 
1/2 fermions plays the main role. Neglect all interac- 
tion terms except the largest on-site one, we introduce 
an effective two-band Hubbbard model on cubic lattices: 
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where a = x, y is the band index, a = j, J, denote spins, 
£ka is energy dispersion for each band. a s is the s-wave 
scattering length which could be tuned by Feshbach res- 
onance technique over a large scope of magnitude. We 
only consider repulsive interaction (a s > 0) in this study. 
Parameter U describes the on-site Coulomb repulsive in- 
teraction between atoms reside in the same orbital, while 
U' describes the repulsive interaction for two atoms in 
different orbitals. J controls the strength of Hund's rule 
coupling between different orbitals. In the present sys- 
tem only transverse part of the Hund's rule coupling sur- 
vives, i.e. the spin flip and pair hopping processes. The 
absence of the density-density interaction for atoms with 
the same spins is due to the short range 5 type pseudo- 
potential. This fact can also be understood as a fully 
canceling between the longitudinal part of the Hund's 
rule coupling and corresponding terms in the inter-orbital 
density-density interactions. All of the interaction pa- 
rameters are expressed as overlaps of maximally localized 
Wannier functions 151. 
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FIG. 2: (a). Quasi-particle weight v. s. interaction strength U. 
U' and J are varied with U by keeping U' = J = 0.299J7. At 
quarter filling (Red hollow dot) the quasi-particle weight satu- 
rate to nonzero value, indicating an absence of the Brinkman- 
Rice transition, (b). Stagger magnetization v.s. interaction 
strength for several different fillings. At half filling (n = 2) the 
system crossover from weak coupling SDW behavior to strong 
coupling spin-1 Neel state. Away from half filling, it's needed 
a finite interaction strength to induce stagger magnetization. 



Both the band dispersions Sk a and the coupling con- 
stants in Eq. [2] are determined by experimentally con- 
trollable parameters V and a s . Due to the anisotropic 
nature of the p-orbitals, the system shows dimension 
reduction behavior: the PDOS has Van Hove singu- 
larity near the band edge resembling of ID character. 
This feature is also included by the tight binding model 
with an anisotropic hopping amplitude along different 
directions . At half filling, the Fermi surface has nest- 
ing property with the nesting vector (jr, n, 7r)(l7|. In con- 
trast to the isotropic case (s-band), here the singular- 
ity in PDOS and Fermi surface nesting occur at separa- 
ble energy scales and lead to remarkable consequence: 
instability towards to different magnetic orders domi- 
nates at different fillings. For the given lattice depth 
V x = V y = 16Er,V z = 36Er, our calculation shows 
the bandwidth W of the two-fold degenerate p-band is 
0.45-Er, while band gaps separated p x , p y -bands with s 
and p z -bands are 6.79E R and 4.10E R (FigrjTb)). Over- 
lap of Wannier functions gives U/Er = 8.63fcia s and 
V = J = 0.299J7, while the interaction scale within s- 
band is 13.1(kLa s )En. Not exceeding the band gaps pose 
0.0755A as the upper limit for a s . Within this restric- 
tion, one could still tune the system from weak to strong 
coupling (U ~ 9W) region by Freshbach resonance tech- 
nique. 

Gutzwiller variation results: For this complex two- 
band Hubbard model we adopt a generalized version of 
the Gutzwiller variational approach [H [Hj] , 
Historically, the variational approach has been success- 
fully used to treat strong correlated systems such as 
normal 3 He [23| and Mott transition Generaliza- 
tion of it is used to investigate physics associate with 
band degeneracy, such as orbital selective Mott transi- 
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tion 2(| 2f| and transition metal ferromagnetism 21 1. 
The variational wave function is |<I>g >= Pg |Q >j where 
|0 > is uncorrelated state which could be fermi liquid, 
spin density wave or superconductivity ground state. 
And Pq = Yii ^i.v m i.Y, mi,r are projection oper- 
ator onto atomic configurations, Aj ; r are corresponding 
variational parameters. Under Gutzwillcr approximation 
(GA) [li| [23] one can evaluate expected values of opera- 
tors over the projected wave function. By minimizing 
the total energy one gets the ground state configura- 
tions, from which one can identify various orders e.g., 
charge density wave order, orbital order, FM and AF or- 
der etc. GA neglects spatial correlations and is only exact 
in infinite dimension, this approximation has also been 
proven to be equivalent to the saddle point of Kotliar- 
Ruckenstein slave boson functional treatment (26J. Since 
it is non-perturbative in nature, the variational approach 
treats the Fermi liquid as well as Mott localized state on 
equal footing and hence gives a coherent description of 
the intermediate coupling strength, connecting the weak 
coupling mean field to the strong coupling perturbation 
results [17] . And the result of the present method is supe- 
rior to Hartree-Fock mean field that often overestimates 
the tendency towards the ordered phases. 

Disregard any long range order, at commensurate fill- 
ing a reasonable large on-site repulsive interaction would 
localize fermions and drive the system into the paramag- 
netic Mott insulator phase. The transition (Brinkman- 
Rice transition) could be described naturally within 
Gutzwillcr variational approach 24] via quasi-particle 
weight Z, which contributes to the coherent peak in the 
spectral function. It could also be interpreted as the 
renormalization factor of kinetic energy, the height of 
Fermi step or inverse of mass enhancement factor. As 
Z approaches to zero, fermions becomes more and more 
localized and finally transforms into the Mott insulator 
phase. However in the present case the large Hund's rule 
coupling J — U' has dramatic effect on the paramag- 
netic Mott transition: there is no Brinkman-Rice transi- 
tion at quarter filling, see Fig[2]Ja). The on-site interac- 
tion Hamiltonian EqQ] has doubly occupied spin-triplet 
as its lowest energy states, which are degenerate with 
the empty and singly occupied states. Charge fluctu- 
ation through this channel is allowed at quarter filling 
thus the system does not become localized when U in- 
creases. At half filling this kind of fluctuation is blocked 
by the large on-site intra-orbital repulsive interaction 
U and the Brinkman-Rice transition manifests itself as 
shown in Fig[2fa). However because of the huge spin 
entropy of the paramagnetic Mott state, it would not 
act as the ground state at zero temperature. When the 
system cools down the degeneracy will be lifted and var- 
ious magnetic/orbital orders develop depending on the 
residual interactions between the spin/orbital degrees of 
freedom. 

Next we discuss the magnetic order of the p-band 
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FIG. 3: Ground state phase diagram. Region enclosed by full 
and dashed line is antiferromagnetic and ferromagnetic phases 
respectively, the remainder part is in paramagnetism phase. 
Contrary to the single band case, there is robust ferromagnetic 
phase for vast range of filling and coupling strength. 



fermions. As the previous section indicated, the band 
structure (Fermi surface nesting and the Van Hove singu- 
larity) favors AF and FM order at different energy scales, 
and the on-site Hund's rule coupling may help stabiliz- 
ing them. Thus, it is possible that these two phases may 
exist even for weak or intermediately coupling strength. 
And it is interesting to observe the coexistence of them 
in the optical lattices with an external harmonic trap. 

Fig. [^b) shows the development of the stagger mag- 
netization with increase of U near half filling (n = 2). 
Due to the Fermi surface nesting, an arbitrary small U 
drives the instability of the two-sublattice antiferromag- 
netism at exact half filing. At weak coupling it follows 
the spin density wave mean field prediction, the stag- 
ger magnetization increases as te~ l / u . Approaching the 
strong coupling limit the stagger magnetization comes to 
its saturation value 2. In this limit, the large repulsive in- 
teraction quenches the charge degree of freedom while the 
on-site Hund's rule coupling locks the local spins to form 
spin-1 moments. Virtual hopping process leads to the 
antiferromagnetic coupling of these local moments and 
the original fermions model reduces to a spin-1 antifer- 
romagnetic Heisenberg model. This high spin AF order 
is more stable against the quantum fluctuations. The 
present method also gives a coherent description of the 
weak to strong crossover region, where the stagger mag- 
netization shows non-monotonic behavior for n = 1.8 as 
shown in Figf^b). The condensation energy of AF or- 
der is roughly W 2 /U, which is 10 times larger than the 
s-band, see FigQJb)(c). Due to the large energy gain of 
the antiferromagnetic order, we anticipate the transition 
temperature to be two orders higher than that of the s- 
band case at the strong coupling limit. At intermediate 
coupling strength the transition temperature may attain 
its maximum value, which is accessible in current cold 
atomic experiments. 

Upon doping, the AF order is destroyed by movement 
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FIG. 4: (a). Side view of spin density profile in a harmonic 
trap with different total number of atoms, two spin state 
are 50 : 50 mixture. The frequency of the external trap is 
Q = 0.005i?R. (b). Top view of various phases. Dark Blue 
region denotes antiferromagntic phases, Black region denotes 
band insulator and shallow blue/red for two kinds of polarized 
fermions (ferromagnetic phase). Only density of atoms on the 
p-band are shown in the figure, s-band spin densities provides 
a homogeneous background near the center of the trap. 



of holes and the system has featureless PM ground state. 
But at even lower fillings where the Van Hove singularity 
plays a role the system has FM ground state. Actually 
for the present depth of optical lattice, the singularity in 
DOS plays the main role and the Hund's rule coupling 
further stabilizes the tendency and enlarges the region of 
FM order in the phase diagram. It is possible to tune the 
optical lattice to change the relative weight of contribu- 
tion of singularity in DOS and the Hund's rule coupling, 
and this will shed light on the long controversial issue of 
the mechanics of ferromagnetism in the transition metal 
oxides. Due to the approximately particle-hole symmetry 
of the DOS, the phase diagram shows a mirror symmetry 
with the n — 2 axis, i.e., there is a FM region of similar- 
shape at higher fillings. 

We now come to the effect of the external confinement. 
It is provided by a combination of the magnetic trap 
and the Gaussian profile of the laser beams. Combing 
with the local spin density approximation(LSDA), our 
Gutzwiller variational approach provides a natural way 
of incorporating spatial inhomogeneities into the treat- 
ment of strong correlations. According to the celebrated 
Hohenberg-Kohn theorem [13], the ground state energy 
of an inhomogeneous system is solely determined by its 
ground state spin density distribution. 
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where E km+mt [ni (r) , n^(r)] includes the kinetic en- 
ergy, the lattice potential and the interaction energy 
of fermions on an optical lattice, which depends on 
the spin density distribution non-locally in general. 



F es *(n T (r),nj.(r)) = ±ft 2 r 2 (n T (r) + nj.(r)) is the trap 
energy. Usually f2 is much smaller compare to the 
characteristic frequency of the optical lattice, thus the 
trap provides a spatially slowly varying chemical po- 
tential. We could expand the functional to the lowest 
order: E m+int [n^(r), n;(r)] « E km+mt {n^(r), n[ {v)), 
where the functional dependence on local spin densities 
could be obtained by previously Gutzwiller variational 
calculation. Physical picture of this strategy (LSDA) is 
to divide the whole system into mesoscopic clusters of 

size Iq = \J jj^. Each of them feels homogeneous exter- 
nal potential and at the same time could be treated as 
in the thermodynamic limit, and Gutzwiller variational 
study gives the dependence of energy and local spin den- 
sities within each cluster. 

We perform optimization of E tot by varying over spin 
density distributions for U = O.GEr (U' and J are de- 
termined correspondingly) and several total number of 
atoms, the results are shown in Fig. |4] Densities of both 
spins decrease from the trap center continuously to zero 
at the trap edge because of the external potential. Due 
to the orbital degeneracy, maximum filling at the trap 
center is 2 for each spin, forming band insulator phases. 
Competition of the external potential and cohesive en- 
ergy of the ordered phase gives fine structures of spin 
density profiles. Regions where spin density around 1 
are more likely to enter into AF phase, as long as the en- 
ergy gain from AF phase exceeds correspondingly energy 
loss from density redistribution. Plateaus formed by AF 
(or PM Mott Insulator) phase was also reported in pre- 
vious Quantum Monte Carlo [28[ and Dynamical Mean 
Field theory [29[ study on single band inhomogeneous 
cold fermionic atom systems. For regions with average 
filling less than 0.5 or larger than 1.5 there are tendency 
towards ferromagnetic order. But due to the conserva- 
tion of total spins, FM region will consist of two do- 
mains of half-shell shape with opposite spin polarization, 
see FigHJa), similar structure is anticipated in [3(| con- 
cerns Nagaoka ferromagnetism. To sum up, p-band cold 
fermionic atoms will form shell structures with an exter- 
nal trap potential, for different radii the system crosses 
the U — n phase diagram and shows an band insulator, 
PM, FM and AF phases Fig. gfb). 

Experimental signatures: Population of higher bands 
have been detected by time of flight (TOF) images Q 31 1 . 
Noise correlation (321 [33| from TOF images may also de- 
tect spin order through spin-spin correlation functions. 
Specifically, AF phase opens a charge transfer gap which 
could be detected by Raman spectroscopy [34j , while the 
doubling of unit cell might be detected by spin selec- 
tive Bragg spectroscopy [35j . Spatial distribution of spin 
densities in harmonic trap reported in this paper could 
be detected by spatially microwave transition and spin- 
changing collisions techniques, which measure the inte- 
grated density profiles along chosen directions [36[ . 
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Summary: Combining the band structure calculations 
and the Gutzwliller variational approach, we perform a 
first principle calculation of the zero temperature phase 
diagram of spin- 1/2 cold fermionic atoms on the two- fold 
degenerate p-band in an optical lattice. We show that 
the system has robust ferromagnetic and antiferromag- 
netic ground state at different fillings. We have traced 
back the physical picture to the single particle feature in- 
cluded anisotropy of orbital orientations and Fermi sur- 
face nesting as well as correlated effects such as Hund' 
rule coupling. We also discuss the inhomogeneous spatial 
distribution induced by an external harmonic potential. 
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